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Using numerical simulations, we show that a novel spontaneous transverse response can appear 
when a longitudinal drive is applied to type-II superconductors with honeycomb pinning arrays in 
a magnetic field near certain filling fractions. This response is generated by dynamical symmetry 
breaking that occurs at fields away from commensurability. We find a coherent strongly amplified 
transverse switching effect when an additional transverse ac drive is applied. The transverse ac drive 
can also be used to control switching in the longitudinal velocity response. We discuss how these 
effects could be used to create new types of devices such as current effect transistors. 

PACS numbers: 74.25.Qt 



There have been extensive studies on superconducting 
systems with patterned pinning arrays which show com- 
mensurability effects when the vortex density matches 
the pinning site density [l], 0, H, 0| ■ Various properties 
of the vortex transport can be controlled by adjusting 
the shape and geometry of the pinning arrays, producing 
vortex channeling [5( , dynamical transitions [a |7| , fluxon 
ratchets [§] , reversible ratchets 0, [13] , and vortex cellu- 
lar automata [HI]. These studies indicate that vortices 
interacting with patterned substrates may lead to a new 
field of fluxtronics or microelectronic devices based on 
the controlled motion of vortices. Further, vortices inter- 
acting with periodic pinning arrays also exhibit a number 
of collective dynamic behaviors which are important to 
the broader field of noncquilibrium physics. 

Few studies have been performed on vortex dynamics 
in honeycomb pinning arrays since it was assumed that 
the vortex behavior in this pinning geometry would be 
similar to that found in square and triangular pinning 
arrays. In this Letter we demonstrate that in fact, hon- 
eycomb pinning arrays produce new types of phenom- 
ena that do not occur in triangular or square pinning 
arrays. In particular, an effective dimerization or higher 
order n-merization of the interstitial vortices in the hon- 
eycomb pinning arrays for certain field ranges results in 
a novel spontaneous transverse response (STR) to a lon- 
gitudinal drive. Normally, an applied current produces 
a perpendicular Lorentz force on superconducting vor- 
tices, which move and generate a voltage drop parallel to 
the current. We show that the honeycomb system can 
also produce a voltage drop perpendicular to the applied 
current. For magnetic fields away from commensuration, 
the vortex flow dynamically organizes to a STR state. 
At the matching fields, the STR arises due to symmetry 
breaking in the ground state. The vortex n-mer states 
are very similar to the colloidal molecular crystal states 
studied for colloids interacting with periodic optical trap 
arrays [12]. We specifically consider vortices in super- 
conductors; however, we expect similar phases to occur 
whenever there is a spontaneous symmetry breaking of 



effective n-mer states for particles on periodic substrates. 

We simulate a two-dimensional system containing N v 
vortices and N p pinning sites with periodic boundary 
conditions in the x and y directions. The vortex den- 
sity is B = N v cf>o/L 2 , where L = 24 A is the system size 
in units of the penetration depth A and cf>o = h/2e is the 
elementary flux quantum. The motion of a single vortex 
is given by the following overdamped equation: 

= F™+Ff + F dc + F ac + Ff. (1) 

Here R^ is the location of vortex i, rj = </>§d/27r£ 2 is 
the damping constant, d is the thickness of the supercon- 
ducting sample, £ is the coherence length, and pjv is the 
normal state resistivity. The pairwise vortex-vortex in- 
teraction force is F™ = Y^iJj /o-^i(-Ry/A)Ry where K\ 
is the modified Bessel function, = |Rj — R,-|, Ry = 
(Rj — Hj)/Rij, and fo = (/>o/(27r^ A 3 ). The pinning force 
F^ arises from parabolic traps of radius r p — 0.3A and 
strength f p = 1.0/o arranged in a honeycomb lattice with 

= - Ef=i f P Rikr- 1 Q((R P - J^ fc )/A)R ife . Here R lk = 
|Rj — Rfc| is the distance between vortex i and pin k, 
R,:fe = (Ri — Rfc)/i?ifc, and 9 is the Heaviside step func- 
tion. The dc longitudinal driving force F dc = F rfc /oR i , 
where we take the longitudinal direction R L to be R = x. 
This force mimics the Lorentz force created by an applied 
current. We also consider the effect of adding a transverse 
ac drive F ac = F ac f sin(wt)R Tr , where w = lCr 7 /oA? 
and R Tr = y is the transverse direction. The thermal 
force Fj arises from Langevin kicks and has the proper- 
ties (*f (t)) = and {FT (t)Fj (t 1 )) = 2r)k B T6 ij 5(t - t 1 ). 
Vortex positions are initialized with simulated annealing. 

In Fig. 1 we illustrate the vortex and pinning site 
positions in a sample with F dc = 0, F ac = 0, and 
T = 0. At B/B^ = 2.0, shown in Fig. l(a,b,c), where 
B^ is the field at which N v — N p , the large interstitial 
space of each plaquette in the honeycomb pinning lat- 
tice traps two interstitial vortices which form an effective 
rigid dimer with a director field that can point in one of 
three degenerate directions, as indicated by the dashed 
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FIG. 1: Vortex positions (black dots) and pinning site lo- 
cations (open circles) in a 20A x 20A section of the sam- 
ple at F dc — 0. (a,b,c): Three different realizations of 
B/B^ — 2.0 with interstitial dimers aligned (a) in the y- 
direction; (b) at +30° to the ^-direction; (c) at —30° to the 
x-direction. The dashed lines indicate the alignment direc- 
tion, (d) BjB^ — 1.77 with no overall dimer alignment. 

lines. Neighboring dimers interact through an effective 
quadrupole interaction which causes the dimers to align 
into one of three degenerate ground states For fill- 
ings 1.5 < B/B^ < 2.0 and 2.0 < B/B^ < 2.5, the dimer 
alignment is disrupted and there is no global symmetry 
breaking in the ground state. Figure 1(d) shows an ex- 
ample at B I Bcj, = 1.77 where the dimer ordering is lost. 

We next apply a dc drive F dc in the longitudinal di- 
rection and measure the vortex velocity v to obtain the 
longitudinal velocity Vl = Yli" v * ' & L an d the trans- 
verse velocity Vr r — J2i * v i ' R- Tr - In Fig- 2(a) we 
plot Vl and Vr r versus F dc for the system in Fig. 1 at 
BjB^ = 2.0. For F dc < 0.14 the system is pinned, while 
for 0.14 < F dc < 0.37 there is a finite longitudinal veloc- 
ity Vl accompanied by a finite transverse response Vr r 
that can be in either the positive (V^ r ) or negative (V^,) 
transverse direction. The trajectories of the moving vor- 
tices for the Vp r and Vf r states are shown in Figs. 3(a) 
and (b), respectively. When F dc > F dc , where F dc is the 
critical drive at which vortices in the pinning sites begin 
to depin, the symmetry breaking is lost. At B/B^ = 2.0, 
F dc = 0.37, and above this drive, |Vrv| drops abruptly to 
zero and Vl rapidly increases. The curve with V^ r was 
obtained by starting from the ground state in Fig. 1(a) 
at T = while the curve with Vf r was generated start- 
ing from the same ground state but applying a finite but 




FIG. 2: (a) Transverse velocities V^ r and longitudinal veloc- 
ity V L versus F dc for the sample in Fig. 1 at B/B4, = 2.0. (b) 
\Vtt\ vs F dc for the same system at B/B^ = 1.58, 1.61, 1.67, 
1.72, 1.77, and 1.81, from bottom to top. (c) V^, ax versus 
B/Bj,. (d) \V T t\/Vl versus B/B^ at fixed F dc = 0.225. (e) 
Time series of V Tr for F dc = 0.3 at B/B^ = 2.0 (top curve) 
and B/B^, = 1.81 (two lower curves), (f) Vrr during consec- 
utive cycles of F dc from to F dc = 0.3 at B/B$ = 2.0 (open 
squares) and B/B^ = 1.81 (filled circles). 



small temperature T = 0.2T m , where T m is the vortex 
lattice melting temperature. If we repeat the finite tem- 
perature simulation with different random seeds, we are 
equally likely to observe V^ r or Vf r . If the initial ground 
state already has a global symmetry preferred positive 
or negative orientation with respect to the transverse di- 
rection, as in Figs. 1(b) and (c), the transverse response 
is in the same direction indicated by the dashed lines in 
Fig. l(b,c). The appearance of a transverse velocity at 
B I Bfj, = 2.0 is not a manifestation of a dynamical sym- 
metry breaking since the dimer alignment symmetry is 
already broken as in Fig. 1. In contrast, for the fields 
1.5 < B/B^ < 2.0 and 2.0 < B/B^ < 2.5 there is no 
global symmetry breaking in the ground states; however, 
a dynamical symmetry breaking occurs when the vortices 
organize into a dynamical phase resembling the states in 
Figs. 3(a,b) after a transient period of time during which 
the vortices move in both the Vp r and Vf r directions. 

In Fig. 2(b) we plot \Vr r \ as a function of F dc at 
B/B^ = 1.58, 1.61, 1.67, 1.72, 1.77, and 1.81 show- 
ing the appearance of the dynamical broken symmetry 
state. In Fig. 2(e) we plot V Tr vs time for F dc = 0.3 at 
BjB^ = 2.0 and B/B^ = 1.81, where the latter curves 
were generated from samples in slightly different initial 
states. At B/B^ = 2.0, the vortices immediately move in 
the broken symmetry direction of the ground state, which 
in this case is V£ r shown in Fig. 1(b). For B/B^ = 1.81, 
the initial motion is symmetric and only develops into a 



3 













; . ; . ; W ; . ; 




(d) X 



FIG. 3: Vortex positions (black dots), pinning site locations 
(open circles), and trajectories (lines) over a constant time 
interval, (a) The V^ r dynamically broken symmetry state 
from Fig. 2(a) at F dc = 0.225. (b) The V^ r state at F dc = 
0.225. (c) A sample with B/B^ = 2.11 at F dc = 0.4 where 
complete symmetry breaking is lost, (d) The V^ r state in a 
sample with B/B^ = 4.5 and F dc = 0.15. 
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FIG. 4: (a) V Tr (t) (thick lower line) and V L (t) (thick upper 
hne), with the analytic values of V ac (t) (thin lower line) and 
Vdc{t) (thin upper line) shown for comparison. Vl and Vdc 
have been shifted up by 0.2 for presentation purposes. Time 
is given in units of 2tv/lj. Here B/B^ = 2.0, F ac = 0.125, and 
F dc = 0.28. (b) Time series of V Tr for a system with B/B^ = 
1.77, F dc = 0.28, and varied ac drive amplitude of F ac = 0.15, 
0.1, 0.075, 0.05, 0.035, and 0.02, from top to bottom on the 
positive response side, (c) The amplification factor A versus 
F ac for samples at F dc = 0.28 with B/B4, = 1.61 (filled 
circles), 2.0 (open squares), and 2.11 (filled triangles), (d) 



|Vty| vs T/T m at BjB^ = 2.0 and F d 



0.3. 



Vrp r or Vrf r broken symmetry state over time. The sys- 
tem may even fluctuate significantly in one direction be- 
fore locking into the other direction. This shows that the 
moving states can have global symmetry breaking even 
when the equilibrium ground state does not. Due to the 
initial symmetric transient moving state, the memory of 
the pinned state configuration is lost at incommensurate 
fields but retained at commensurate fields when F dc is 
repeatedly cycled from zero into the moving broken sym- 
metry state and back again. In Fig. 2(f), the plot of Vtt 
at B/B$ = 2.0 and 1.81 during consecutive cycles of F dc 
shows that the sign of the symmetry breaking switches 
randomly between cycles at the incommensurate field. 

In Fig. 2(c) we plot \j^ ax , the maximum value of 
|Vt>| obtained when varying F dc , versus B/B^. The 
overall maximum amount of transverse motion occurs at 
B/B<p = 2.0. The fraction of the velocity that is in the 
longitudinal direction |Vrr|/Vt is shown in Fig. 2(d) for 
fixed F dc — 0.225 as a function of B/B^. For complete 
symmetry breaking, the vortices move at ±30° from the 
longitudinal axis so we expect |V^y|/Vx, = tan(30°) = 
0.577, as seen for 1.5 < B/B^ < 2.0 in Fig. 2(d). For 
B/Bcj, > 2.0, both V^ r ax and \V Tr \/V L drop rapidly. 

The same type of symmetry breaking flow occurs at 
higher commensurate states such as B/B^ = 4.5 where 



elongated trimers can form, as illustrated in Fig. 3(d). 
This is accompanied by spontaneous dynamical symme- 
try breaking over a range of fields. We have also con- 
sidered the role of temperature. Previous work (l3| indi- 
cated that the interstitial vortices melt at a well defined 
temperature T m ; using parameters appropriate for Nb 
crystals gives T m w 8.5 K. We find that our results are 
robust for all T < T m . In Fig. 4(d) we plot \V Tr \ vs T/T m 
at B/B& = 2.0 and F dc = 0.3, showing that Vr r drops 
sharply with increasing T just below T/T m = 1.0. 

Since the dynamical symmetry-broken states are 
bistable, it is natural to ask whether it is possible to 
induce a switching behavior between the two states. To 
address this, we consider the effect of adding a trans- 
verse oscillating drive F ac to a system with B/B^ = 2.0 
moving under a fixed longitudinal drive F dc — 0.28. 
Experiments have already been conducted in supercon- 
ductors with periodic pinning arrays where two orthog- 
onal driving currents were simultaneously applied and 
the transverse and longitudinal responses were simulta- 
neously measured 0. In Fig. 4(a) we plot the time 
series Vrv(t) and Vl(£) in a system with F ac = 0.125. 
Here the ac drive induces a periodic switching between 
the positive and negative transverse response states V^> 
and Vf r . We find that this abrupt switching effect per- 
sists for 1.5 < B/B^ < 2.0. Similar switching appears 
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over a wide range of ac driving amplitudes, as illustrated 
in Fig. 4(b) where we plot V Tr (t) at B/B^ = 1.77 and 
pdc _ o.28 for different values of F ac . We observe a novel 
response in the longitudinal direction in which the value 
of Vt(t) varies by up to 50% during each ac drive period. 
There is also a switching effect in Vi(t) which accompa- 
nies each switch in Vr r {t)- This phenomenon is reminis- 
cent of the current effect transistor found in transversely 
driven charge density waves where a transverse force can 
be used to control the longitudinal response [IH . 

There is an amplification of the transverse ac response 
Vr r {t) over the value expected based only on the magni- 
tude of the ac input driving signal. To illustrate this, in 
Fig. 4(a) we plot V ac (t) = N m F ac (t) /rj, the transverse ac 
velocity produced when the N m — N v — N p interstitial 
vortices move only in response to the ac driving force, 
along with Vd c = N m F dc /-q, the longitudinal velocity of 
the interstitial vortices under only the dc drive. Fig. 4(a) 
indicates that |Vr r (i)l > \Vac(t)\- The relative amplifi- 
cation VTr{t)/V ac {t) can be increased by lowering F ac ; 
however, for a finite F dc there is a threshold value F^ c 
below which the system no longer switches between V^ r 
and Vrf r . For example, in Fig. 4(b), F ac — 0.02 is below 
the threshold F£ c , and thus the response stays locked in 
the Vjt direction. We quantify the amplification of the 
transverse ac response using A = Vr r iJ / (N m F ac ) , where 
Vtt is the amplitude of the transverse response, such that 
for A — 1.0, there is no amplification. The plot of A vs 
F ac in Fig. 4(c) for B/B^ = 1.61, 2.0, and 2.11 shows 
that for a fixed dc drive, the amount of ac amplification 
that occurs depends on both F ac — F c ac and F dc — F dc . 
The closer the ac and dc drives are to the respective crit- 
ical thresholds, the larger the amplification. Since F dc is 
maximized at B/B^ = 2.0, for fixed F dc — 0.28 a larger 
amplification can be obtained at B/B^ ^ 2.0 than at 
B/Bcf, = 2.0, as shown in Fig. 4(c). For the parameters 
used in Fig. 4, the frequencies for a typical superconduc- 
tor are between 10 and 30 Hz. If to is increased, the sharp 
switching response can still be achieved by increasing F ac 
until io ~ 20 to 60 KHz, at which point F c ac crosses above 
the depinning threshold for the pinned vortices and the 
dynamical symmetry breaking is completely lost. 

Many of the switching features, including bistability in 
the transverse response and the abrupt switching effect 
induced by the ac drive, are analogous to semiconduc- 
tors and may be useful for creating logic devices. The 
strong coupling between the longitudinal and transverse 
responses indicates that a transverse ac drive could be 
used to control longitudinal switching, while the ampli- 
fication effect implies that very small ac inputs can in- 
duce large response changes in certain regimes. Our re- 
sults should apply to any system of repulsively interact- 
ing particles confined by honeycomb pinning arrays. We 
have observed a similar dynamical symmetry breaking 
flow for colloids on triangular substrates [l6j |: however, 
the dynamical symmetry breaking for the vortex system 



presented here is much more robust and occurs for a con- 
siderably wider range of parameters and particle fillings. 

In summary, we have demonstrated that vortices in 
superconductors with honeycomb pinning arrays can ex- 
hibit a novel transverse response when a longitudinal 
drive is applied for certain ranges of fields where the 
vortices in the large interstitial regions of the pinning 
array form effective dimer or higher order n-mer states. 
At commensurate fields, the symmetry is broken in the 
ground state, while at incommensurate fields, there is a 
dynamical symmetry breaking. If a transverse ac drive is 
added to the system, a pronounced transverse switching 
response occurs, as well as an amplification of the trans- 
verse ac signal. There is a threshold ac drive required to 
induce the switching which becomes very small when the 
longitudinal drive is close to the value above which the 
transverse response disappears. The transverse ac drive 
can be used to modulate the longitudinal response of the 
system as well. We discuss how this effect may be useful 
for creating new types of fluxon based devices. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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